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ABSTRACT 

We report the discovery of two radio-loud quasars with redshifts greater than four; 
GB1428+4217with z = 4.72 and GB1713+2148 with z = 4.01. This doubles the number 
of published radio-selected quasars with z > 4, bringing the total to 4. GB1428+4217 
is the third most distant quasar known and the highest redshift radio and X-ray source 
currently known. It has a radio flux density at 5GHz of 259 ± 31 mJy and an optical 
magnitude of R ~ 20.9. The rest frame absolute UV magnitude, M„(1450A), is —26.7 
similar to that of the archetypal radio-selected quasar 3C273 (z — 0.158; M ly (1450A) = 
-26.4). GB1428+4217 has a tentative detection in ROSAT PSPC observations, which 
has been confirmed by more recent ROSAT observations, described in a companion 
paper by Fabian et al. Both quasars were discovered during the CCD imaging phase of 
an investigation into the evolution of the space density of radio-loud quasars at high 
redshift. Combined with our earlier survey results these objects give a lower limit on the 
space density of quasars with radio power P5GHZ > 5.8 x 10 26 WHz _1 sr _1 between z — 4 
and z = 5 of 1.4±0.9 x 10~ 10 Mpc- 3 . This can be compared to 2.9±0.2 x 10~ 10 Mpc- 3 at 
z = 2 from Dunlop & Peacock (1990) for flat-spectrum sources of the same luminosity. 
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1 INTRODUCTION 



In two recent papers we have described results from the first 
phase of a survey for high-redshift, radio- loud quasars (Hook 
et al 1995, 1996). Our approach involves the optical identi- 
fication of flat spectrum radio sources and the spectroscopic 
follow-up of the red stellar identifications. This approach 
exploits the observation that quasars at high redshift have 
redder optical colours than their low-redshift counterparts 
due to absorption by intervening HI (see figure 1 in Hook et 
al. 1995). 

Our aim is to study the evolution of the quasar pop- 
ulation at high redshift, using a well-defined, statistically 
complete sample. Whilst radio-loud quasars are only a small 
subset of the quasar population, the selection of radio-loud 
quasars is less prone to selection effects than optical sam- 
ples, since the spectral energy distribution of radio sources 
is smooth and radio emission is unaffected by either intrinsic 
or extrinsic absorption due to dust. See Fall & Pei (1995) 
for a discussion of how dust within intervening galaxies may 
affect the observed evolution in optically selected samples of 
quasars. Our work complements the numerous investigations 
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into the redshift evolution of the space density of radio quiet 
quasars (eg. Hall et al., 1996; Hawkins & Veron, 1996; Irwin, 
McMahon & Hazard, 1991; Kennefick et al., 1996; Schmidt, 
Schneider & Gunn, 1995; Warren, Hewett & Osmer, 1995). 

Our earlier work was restricted to sources that were 
detected on APM scans of POSS-I plates and hence was 
incomplete for optically fainter quasars. Since it is our aim to 
determine the evolution of the space density of high-redshift 
quasars independent of optical selection effects, we have now 
begun CCD imaging of sources that were not detected on 
the POSS-I plates. 

Throughout this paper we have assumed cosmological 
constants of Ho = 50kms _1 Mpc -1 and go = 0.5. 



2 THE RADIO SAMPLE AND THE OPTICAL 
IDENTIFICATION PROCEDURE 

The parent radio sample used in this study consists of ~1600 
flat-spectrum (af^QH Z > —0.5, S cx u a ) radio sources with 
S5GHZ > 0.2Jy studied with the VLA by Patnaik et al. 
(1992). Initial optical identification of the radio sources was 
carried out using the Automated Plate Measurement (APM) 
Facility at Cambridge, U.K. Further details of this program 
and the procedures can be found in Hook et al (1995,1996). 
For the second phase of the survey, we have concen- 
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trated on the radio sources with no APM/POSS-I iden- 
tification within a search radius of 3.0". The area stud- 
ied lies in the region with 9 h 50' < a(B1950) < 17 h 45', 
20° < <5(B1950) < 75°, (|6| > 20° for h to I2 h and \b\ > 30° 
for 12 to 24 h ), an area of 1.27sr containing 526 sources. 
There were 82 blank fields with Ssghz > 0.2Jy in this re- 
gion, after our earlier work. Of these, 16 had redshifts from 
other sources prior to the start of this project (NEr|j Ver- 
meulen et al 1996, Vermeulen, private communication). 

An unbiased sub-sample of 38 of the remaining 66 
sources were imaged during one week (29 May to 4 June 
1995) at the Prime Focus of 2.5m Isaac Newton Telescope, 
La Palma. We used the 1024 x 1024 thinned TEK CCD with 
a projected pixel size of 0.59". Mould B and R filters (The 
'Kitt-Peak' set at La Palma) were used with 300s exposures 
in R and 600s exposures in B, and typically reaching limiting 
magnitudes of R = 23mag, B = 24mag. 

Photometric calibration was carried out using observa- 
tions of the photometric standard star fields from Landolt 
(1992). The data were reduced using standard procedures 
with the IRAF0 software environment including photome- 
try using the DAOPHOT package (Stetson, 1990). 

After this INT imaging, 10 sources are still blank in the 
R band, i.e. they had no optical counterpart within a radius 
of 3" (K band imaging of these sources is being obtained). 
A colour magnitude diagram for the 28 objects detected in 
R is shown in Figure |l[ 

Thus after the INT imaging the fraction of sources ob- 
served is 38/66 = 0.58. To estimate in a fair manner the 
effective area observed, we include the same fraction (9 out 
of 16) of the sources with previously known redshifts. The 
effective area is then 0.58 x 1.27sr = 0.73sr. Scaling from 
the original 526 radio sources in 1.27sr, there are 302 radio 
sources in this sub-area of which all but 10 (3%) are now 
identified. 



3 SPECTROSCOPIC OBSERVATIONS AND 
RESULTS 

We have selected for spectroscopic followup a complete sam- 
ple of the 12 stellar identifications with B — R > 1.0. 
Redshifts became available for 3 of the 12 selected ob- 
jects prior to the spectroscopic phase of this study (Ver- 
meulen, private communication) including the z — 3.82 
quasar GB 1239+3736 (Vermeulen et al 1996). 7 of the other 
9 were observed at either the Shane 3m telescope at Lick ob- 
servatory or the 4.2m William Herschel Telescope (WHT), 
La Palma. The Lick spectra were obtained using the KAST 
double spectrograph at lowest resolution covering the range 
5000 - 10500A. The WHT spectra were obtained using the 



t The NASA/IPAC Extragalactic Database (NED) is operated 
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t IRAF is distributed by the National Optical Astronomy Obser- 
vatories, which is operated by the Association of Universities for 
Research in Astronomy, Inc. (AURA) under cooperative agree- 
ment with the National Science Foundation. 
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Figure 1. Colour-magnitude diagram from INT CCD identifica- 
tions for the SsGHz > 0.2Jy flat-spectrum radio sample. Filled 
circles show unresolved sources, open circles represent galaxies. 
The squares show spectroscopically-confirmed 2 > 3 quasars. The 
dashed line indicates the boundary used for the spectroscopic sub- 
sample. The filled triangles show the 2 objects of in that sample 
for which a spectrum has yet to be obtained (see text). 



red arm of the ISIS double spectrograph. A 158 1/mm grat- 
ing and a thinned TEK CCD with 1024x1024 24^m pixels 
was used, giving a useful wavelength range of 5970-8940A 
at 2.9A per pixel. 

Of the 7 objects observed, 1 still has an inconclusive 
spectrum. Another (GB1712+493) has a recently-published 
redshift of 1.552 (Falco, Kochaneck & Munoz 1997), con- 
sistent with the single emission line seen in our spectrum 
being identified with Mgll. The spectrum of one other ob- 
ject also shows a single broad emission line (probably also 
Mgll) which we can rule out as being Ly-a at z > 4 by the 
lack of other strong emission lines or Ly-a forest absorp- 
tion. Here we give details for two objects, GB1428+4217 
and GB1713+2148, which were identified as high-redshift 
quasars. The INT B and R band images of the 2 new z > 4 
quasars are shown in Figure |^. 

GB1428+4217 (z=4.715±0.010) GB1428+4217 was 
first observed on 1996 16th July on the 4.3m WHT and 
then again on 1996 19th July at a redder wavelength grat- 
ing setting in order to detect the CIV emission line. Broad 
emission lines of Lya, NV, SilV/OIV] and CIV are visible 
in the spectrum (Fig. ^). The observed centroids and and 
equivalent widths are given in Table 1. 

The Lya line is visibly asymmetric due to absorption 
and it is therefore impractical to measure the line centroid. 
If one assumes that the absorption is due to intervening HI, 
the blue edge of the line provides a lower limit to the redshift. 
The peak of the blue edge occurs at 6974A. The correction 
for instrumental resolution effects is ~ 3A, so assuming a 
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Figure 2. B and R band images, obtained at INT, of the two z > 4 quasars. The upper panels show GB1428+4217 and the lower panels 
show GB1713+2148. The central cross is 1 arcmin across. North is at the top, West to the right. The offsets (An, AS) in arcseconds 
from the stars marked A and B to the quasars are as follows (positive Aa represents a move East from the offset star to the quasar, and 
positive AS represents a move North): GB1428+4217: from A to quasar (-13.6, 22.3), from B to quasar (80.7 -12.7); GB1713+2148: 
from A to quasar (30.2, -19.9) from B to quasar (-46.0, -14.6). 



value of 697lA for the Ly-a edge gives a redshift of 4.734. It 
is quite common to observe velocity shifts between the emis- 
sion lines in quasars and this redshift difference corresponds 
to a blue shift of 1420 km/sec of CIV with respect to Ly-a. 
This is significantly higher that the mean shift of 138 km/sec 
derived by Tytler and Fan (1992) in their analysis of a sam- 
ple of 160 quasars. We take the redshift of GB1428+4217 to 
be the mean of the three estimates. 

On the night of 19th July, both narrow and wide 
slit (7 arc second) spectrophotometric observations of 
GB1428+4217 were obtained. The flux standard G138- 
31 (1625+0919; Filippenko & Greenstein, 1984) was ob- 
served both prior to and immediately after the observa- 
tions of GB1428+4217. These observations give a contin- 
uum flux at a rest wavelength of 1500A (A o f, s =8550A) of 
6.3x 10 _28 ergs cm _2 s _1 Hz _1 . The corresponding continuum 
apparent magnitude and absolute magnitude on the AB sys- 
tem (Oke & Gunn 1983) are 19.4 and -26.7 respectively. For 



comparison, from HST observations of the archetypal radio- 
selected quasar 3C273 (z = 0.158) (Bahcall et al 1991) we 
find the continuum luminosity to be M„(1450) = —26.4. 

GB1713+2148 (z=4.011±0.005) The discovery spec- 
trum of GB1713+2148 was obtained at Lick-3m on UT 1996 
May 23. It was reobserved at the WHT on 1996 16th July for 
1800s to obtain a higher signal-to-noise spectrum, shown in 
Figure |^. Broad emission lines of Lya, NV, SilV/OIV] and 
CIV are visible in the spectrum. Again the Lya line has 
visible absorption. 



4 DISCUSSION 

GB1428+4217 is the most distant radio-loud object known 
and the 3rd highest redshift quasar known (behind 
PC1247+3406 at z = 4.897 and PC1158+4635 at z = 4.73, 
Schneider, Schmidt & Gunn 1991b, 1989). 
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Table 1. Properties of the new z > 4 quasars. Redshifts were measured as described in the text, assuming rest-frame wavelengths for 
the emission lines as follows: Lya:1215.7, CIV: 1549.1, NV:1240.1, SilV/OIV] : 1397.8 (see Tytler & Fan 1992). The radio fluxes arc from 
Gregory and Condon (1991). The R and B magnitudes are accurate to 0.07mag (rms) or better. *=Obscrved equivalent width measured 
for the Lyce-NV combination, p= wavelength measured from the peak rather than the centroid of line. 
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Figure 3. Optical spectra of the two z > 4 quasars. 



GB1428+4217 has a tentative detection in the ROSAT 
All Sky Survey (Brinkman et al. 1997) and an archival PSPC 
pointing (WGACAT, White Giommi & Angelini, 1994). The 
extremely high apparent X-ray luminosity of this object and 
the possibility that the object is beamed are discussed in 
another Letter (Fabian et al 1997). In Figure 4 we compare 
the equivalent widths of Ly-a+NV for the four known z > 4 
radio-selected quasars with other samples of z > 4 quasars. 
GB1428+4217 is notable in that it is amongst the lowest 
equivalent width objects. This is consistent with a scenario 
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Figure 4. Rest frame Ly-a+NV equivalent width distribution 
for z > 4 quasars. The solid line is for 23 colour selected APM 
BRI quasars from Storrie-Lombardi ct al. (1996). The dotted 
line is for the 10 grism-selected quasars from Schneider, Schmidt 
& Gunn (1991a). The solid bins are for radio selected z > 4 
quasars including both those presented here and GB1508+5714 
and PKS1251-407 (Shaver et al 1996a). 



in which the continuum region is beamed. In addition, in 
contrast to the other three radio loud z > 4 quasars, the 
emission lines are relatively broad. 

After our INT and WHT observations 2 of the 38 POSS- 
APM blank fields are quasars with z > 4 and 23 are ruled 
out as being at z > 4 quasars (either based on spectro- 
scopic observations, or deduced from their blue B — R colours 
or extended appearance on the CCD images). 10 are still 
unidentified and a further 3 of the CCD identifications still 
need conclusive spectra (1 of these has an inconclusive WHT 
spectrum and 2 were not attempted). Clearly, it is impor- 
tant to determine the redshift status of these remaining 13 
sources. Thus in an effective sub-area of 0.73sr containing 
302 sources, 97% (all but 10) are optically identified and 
96% (all but 13) can be classified as being a z > 4 quasar 
or otherwise. 
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Therefore we can derive a lower limit to the space den- 
sity at z > 4 based on our sample, independently of any as- 
sumptions about the form of evolution. In an area of 0.73sr, 
two objects with z > 4 were found among the CCD iden- 
tifications. In our previous work we have found one z > 4 
quasar among POSS-I identifications in an area of 3.66sr 
(Hook et al 1995) which scales to 0.2 objects in 0.73sr. Thus 
a total of 2.2 ±1.4 objects were effectively discovered in 
an area of 0.73sr (where the uncertainties are from Poisson 
statistics) . 

At z = 5 the flux limit of our sample, Ssghz = 
0.2Jy, corresponds to a limiting radio power of 5.8 x 
10 26 WHz _1 sr , assuming a radio spectral index a = 0.0. 
For objects brighter than this limit we derive a space den- 
sity between z = 4 and z = 5 of 1.4 ± 0.9 x 10" 10 Mpc" 3 . If 
any of the 13 objects above are found to lie at z > 4, this 
space density estimate would be increased, hence it should 
be considered a lower limit. 

This can be compared to the value at z — 2 from Dun- 
lop & Peacock (1990) for flat-spectrum sources of the same 
luminosity. By integrating their models 1 — 5 over radio pow- 
ers greater than 5.8 x 10 26 WHz -1 sr -1 , we derive 2.9 ±0.2 x 
10" 10 Mpc" 3 . The study of Shaver et al (1996b) deals with 
somewhat stronger sources (P > 1.1 x lO^WHz-^r -1 ) and 
since they show space densities normalised to z ~ 3, is not 
possible to compare results directly. 
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